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Zero Temperature String Breaking with Staggered Quarks * 
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The separation of a heavy quark and antiquark pair leads to the formation of a tube of flux, or "string" , which 
should break in the presence of light quark-antiquark pairs. This expected zero-temperature phenomenon has 
proven elusive in simulations of lattice QCD. In an extension of work reported last year we present clear evidence 
for string breaking in QCD with two flavors of dynamical staggered sea quarks and apply our results to a simple 
three-state mixing model for string breaking. We find that mixing is weak and falls to zero at level crossing. 



1. INTRODUCTION 

The heavy quark-antiquark potential is known 
quite accurately in quenched simulations Q. It is 
traditionally measured with the Wilson-loop ob- 
servable, proportional to exp[— V(R)T] at large 
T . In the presence of dynamical quarks it is ex- 
pected that the potential levels off with increasing 
R, signaling string breaking. However, in QCD 
simulations string breaking has proven to be very 
difficult to detect with this observable, even out 
to R w 2 fm§,§. 

The reason string breaking has not been seen 
using the traditional Wilson-loop observable has 
been apparent for some time ||-q]- The Wil- 
son loop can be regarded as a hadron correlator 
with a source and sink state "F" consisting of a 
fixed heavy quark-antiquark pair and an associ- 
ated flux tube. The correct lowest energy con- 
tribution to the Wilson-loop correlator at large 
R should be a state "M" consisting of two iso- 
lated heavy-light mesons. However, such a state 
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with an extra light dynamical quark pair has poor 
overlap with the flux-tube state, so it is presum- 
ably revealed only after evolution to very large T. 
To hasten the emergence of the true ground state, 
it is necessary to enlarge the space of sources to 
include both F and at least one M state. 

String breaking has been demonstrated in the 
strong-coupling approximation to QCD ||f7| , and 
a variety of QCD-like theories prO|, and in 
SU(3) with dynamical quarks at nonzero temper- 
ature^]. Last year, we reported a preliminary 
low-statistics result for two flavors of staggered 
quarks E3, and, earlier this year, Pennanen and 
Michael announced evidence for string breaking 
at zero temperature using Wilson-clover quarks 
and a novel technique for variance reduction in 
computing the light quark propagator |ll| . This 
year we report results with higher statistics u& . 



2. METHODOLOGY 



We work with 198 configurations of size 20 3 x 
24, generated with the one-plaquette gauge action 



Figure 1. The static- light meson- ant imeson pair 
contribution to the full QCD propagator. The 
wiggly lines denote the light quark propagator. 
Shown are the 'direct' and 'exchange' terms re- 
spectively. 
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Figure 2. The string-meson correlation matrix el- 
ement Gfm (and its hermitian conjugate Gmf)- 
The wiggly line again denotes the light quark 
propagator. 

at Q/g 2 = 5.415 in the presence of two flavors of 
conventional dynamical staggered quarks of mass 
am = 0.0125. The lattice spacing (via the Som- 
mer parameter ]l7|]) is approximately 0.163 fm, 
and /nip — 0.358, giving a comfortable box 
size and a relatively light quark. 

Our flux tube "F" source and sink operator is 
constructed from the product of smeared spatial 
links, for the most part taken along one of the lat- 
tice axes, using 10 APE smearing iterations jl8| , 
combining the direct link with a factor 1 — a (in 
our case, a = 0.294) and six staples with factor 
a/6 with SU(3) projection after each iteration. 
The correlator of this operator Gff(R,T) is the 
familiar Wilson loop with smeared space-like seg- 
ments, and point-like static quark and antiquark 
lines. Our two- meson source and sink operator 
"M" is the direct product of static-light meson 



Figure 3. Static-light propagator with a nonoscil- 
lating S- wave and oscillating P-wave component. 
The solid line connects the best fit values. 

and antimeson operators .[] In constructing the 
static-light meson we use a light-quark wave func- 
tion with weight 2 at zero separation and 1 on 
each of the six on-axis second neighbors. Thus 
we measure the additional correlation matrix ele- 
ments Gmm(R,T), Gmf(R,T) and Gfm(R,T), 
diagramed in Figs, [j] and g. These observables 
are computed for all translational and cubic ro- 
tational displacements. The required all-to-all 
propagator is estimated using a random source 
technique with 128 random sources per configu- 
ration. 

3. RESULTS 

The static-light meson correlator in Fig. || has 
contributions from both a nonoscillating S*-wave 
(B-meson-likc) light-quark orbital and oscillat- 
ing P-wave (B* -meson-like) orbital with aEs = 



1 In the continuum limit the meson pairs created by this 
operator are in a combination of light quark singlet and 
nonsinglet states, and only the singlet part should mix 
with the flux tube. On our coarse lattice, where only one 
meson-meson state is important, this affects the overlap 
Zm-, but not the mixing parameters. In this analysis we 
have not attempted to reconcile the restriction to two fla- 
vors in internal quark loops with the flavor counting of the 
external states, which would affect the relative weights of 
the "direct" and "exchange" contributions to Gmm- With 
the accuracy we have achieved so far, this makes no de- 
tectable difference. 
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0.7884(12) and aE P = 1.022(6). Thus our two- 
meson correlator includes nonoscillating combina- 
tions SS and PP and an oscillating SP. If mix- 
ing is weak, our two-channel correlators should 
also include a flux-tube level with an approxi- 
mately Coulomb-plus-linear behavior in R. Our 
fit ansatz for the two-channel correlator is 



N 



G AB (R,T) ='£z% i (R)ZB i (R)l* i (R)] 



T+l 



(1) 



where A, B refer to the flux tube F or meson- 
meson M states and Xi(R) is a real (positive or 
negative) eigenvalue of the transfer matrix. The 
channel energy is — log(|Aj(-R)|). 

For simplicity we restrict the analysis to three 
spectral components 



\ X {R) = , 
A 2 (i?) = 

HR) = (-) 



A 2 (i?) = -e-W 



(2) 



R+l e -V 3 (R) 



corresponding in the unmixed language, respec- 
tively, to the SS, SP, and flux-tube components. 
The peculiar phase for the flux tube component is 
required by discrete lattice symmetry and can be 
obtained simply by interpreting the Wilson loop 
as a heavy staggered quark loop. 

Results for the three levels are plotted in Fig. |[ 
We see clear evidence for string breaking at the 
first level crossing (about 1 fm), but with our 
statistics, we see no evidence for rounding asso- 
ciated with avoided level crossing. Thus string- 
breaking pair-creation and annihilation transi- 
tions are weak in QCD. 

To explore mixing further, we introduce a mix- 
ing model similar to that of Drummond, Penna- 
nen and Michael |6|]il|. 

(3) 



G(R,T) = Z°(R)T{R) T+1 Z a (R) , 
where the transfer matrix is given by 
A°(i?) 



T(R) = 





x 



X° 2 (R) 
V 



x 



A§(i?) 



(4) 



and the mixing coefficients are x and y. This 
model matches our fit ansatz reasonably well, al- 
lowing us to extract values for the mixing pa- 
rameters, the magnitudes of which are plotted in 



rp ffl 

— i 



9 -S ! 



10 



R 



Figure 4. Heavy quark potential and first two 
excited states vs separation R. The dashed and 
solid lines give the asymptotic values 2aEs and 
Ep + Eg. Jackknife errors are shown. 
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Figure 5. Absolute value of the mixing parame- 
ter x vs separation R. Odd and even series are 
distinguished. 
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Figure 6. Absolute value of the mixing parame- 
ter y vs separation R. Odd and even series are 
distinguished. 

Figs. H and ^. We see that mixing is weak and is 
small in the vicinity of level crossing. 

4. DISCUSSION 

To see string breaking in the heavy-quark po- 
tential, it is helpful to introduce explicit meson- 
meson channels. Mixing among the channels is 
found to be very weak. String breaking occurs 
at the first level crossing, i.e., with our quark 
mass, about one fermi. These results justify a 
two-channel model of excited quarkonium decay, 
in which a closed channel populated by bound 
states couples weakly to an open channel jl9| . 

This work is supported by the US National Sci- 
ence Foundation and Dcpartmet of Energy and 
used computer resources at Indiana University, 
the San Diego Supercomputer Center (NPACI), 
and the University of Utah (CHPC). 
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